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In order to search for the magnetic ground state of surface nanostructures we extended first 
principles adiabatic spin dynamics to the case of fully relativistic electron scattering. Our method 
relies on a constrained density functional theory whereby the evolution of the orientations of the 
spin-moments results from a semi-classical Landau-Lifshitz equation. This approach is applied to a 
study of the ground state of a finite Co chain placed along a step edge of a Pt(lll) surface. As far as 
the ground state spin orientation is concerned we obtain excellent agreement with the experiment. 
Furthermore we observe noncoUinearity of the atom-resolved spin and orbital moments. In terms 
of magnetic force theorem calculations we also demonstrate how a reduction of symmetry leads to 
the existence of canted magnetic states. 
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Stimulated by the need for ever higher density record- 
ing media, atomic scale magnetic devices are presently 
at the very focus of experimental and theoretical re- 
search (see, e.g., the "viewpoint" drawn by Kiibler 0]). 
Without doubt, understanding and design of the rel- 
evant physical properties - magnetic moments, mag- 
netic anisotropy energies, thermal stability, switching - 
of atomic scaled magnets demand detailed knowledge of 
their electronic and magnetic structure. Several studies 
based on the Hubbard model S Q or on density func- 
tional theory (DFT) HHHSlI have attempted to ex- 
plore the, mostly, noncoUinear spin ground state of free 
and supported metallic clusters. While noncoUinearity 
due to frustrated exchange interactions can be described 
in terms of a nonrelativistic theory 0, H, IE 1 general 
spin structures are subject to an interplay between the 
exchange and the spin-orbit interaction and, therefore, 
can only be studied within the framework of a relativis- 
tic electron theory Q . 

An efficient tool to search for the equilibrium spin ar- 
rangement of a spontaneously magnetized system is to 
trace the time evolution of the spin moments until a sta- 
tionary state is achieved. The foundations of the so- 
called first principles adiabatic spin-dynamics (SD) for 
itinerant electron systems were laid by Antropov et al. 
in quite a general context. In short, for systems with 
well-defined local (atomic) moments, the evolution of the 
time dependent orientational configuration, {ei(t)}, is 
described by a microscopic, quasi-classical equation of 
motion. 
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where '^1^^ is an effective magnetic field averaged over 
the cell z, 7 is the gyromagnetic ratio and A is a damp- 



ing (Gilbert) parameter. In this equation, the instanta- 
neous orientational state is evaluated in terms of a self- 
consistent calculation within DFT. This formalism was 
further developed by Stocks et al. llOl llll employing a 
constrained density functional theory Here a local 
(transverse) constraining field, B™" - that can be deter- 
mined selfconsistently - ensures the stability within DFT 
of the nonequilibrium orientational state demanded by 
the equation of motion. Clearly, the internal effective 
field that rotates the spins in the absence of a constraint 
and, therefore, has to be used in Eq. (Q) is just the op- 
posite of the constraining field ^3 ■ By merging with the 
locally selfconsistent multiple scattering (LSMS) method 
SD has been applied to study the coniplex magnetic or- 
derings in bulk metals and alloys [lOl IllL Il3| . 

In surface nanostructures, spin-orbit coupling - its im- 
portance magnified by the reduced, surface symmetry - 
obviously plays a key role in determining magnetic orien- 
tations. In order to deal with exchange splitting and rela- 
tivistic scattering on an equal theoretical footing, we use 
the above first principles constrained DFT-SD scheme in 
conjunction with the Kohn-Sham-Dirac equation. 
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In this equation, a. and (3 are the usual Dirac matri- 
ces, (T are the Pauli matrices, V{y) includes the Hartree 
and exchange-correlation potentials, while, within the lo- 
cal spin density approximation (LSD A), B^'^(r) is an ex- 
change field interacting only with the spin of the electron. 
It should be noted that Eqs. and ^ display the very 
basics of a relativistic spin-only dynamics, inasmuch no 
attempt is made to trace, distinctly, the time evolution of 
the orbital moments. Furthermore, although implied by 
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use of Eq. ^ with Eq. we did not include a torque 
related to the coupling of the spin- and orbital degrees 
of freedom. As pointed out by Antropov et al. 9] this 
approach is only applicable when the deviation between 
the orientations of the spin and orbital moments is small. 

In order to apply SD to nanostructures of finite size we 
have merged the above scheme with the multiple scat- 
tering theory (MST) Green's function embedded clus- 
ter method developed by Lazarovits et al. 0|. In 
this, a self-consistent calculation is first carried out for 
the surface system in terms of the relativistic Screened 
Korringa-Kohn-Rostoker method 0|. The nanostruc- 
ture is subsequently inserted into this host using the con- 
ventional MST embedding formula, 

^(6) ^ T-(6) [I - (t^ie)-' t%e)-') T-(6)] . (3) 

In this equation, T^(t'') and T^(t^) are site-angular 
momentum matrices of the scattering path operators 
(single-site t-operators) of the host surface system and 
the cluster, respectively, and e is the energy. By solv- 
ing this equation together with the corresponding Pois- 
son equation - with appropriate boundary conditions - a 
selfconsistent calculation for the selected cluster can be 
performed that takes full account of the environment. 

Recently, Gambardella et al. 0| reported the results 
of experiments on well characterized finite linear chains 
of Co atoms located at a step edge of a Pt(lll) sur- 
face terrace. At 45 K the formation of ferromagnetic 
spin-blocks of about 15 atoms was observed with an 
easy magnetization axis normal to the chain and point- 
ing along a direction of 43° towards the step edge. Previ- 
ous DFT based studies related to this experiment calcu- 
lated spin and orbital moments |l7l llSi Il9l| and magnetic 
anisotropy energies of finite|0| and infinite chains on 
Pt(lll) surfaces using a magnetic force theorem (MFT) 
approach "2^. In what follows, we report the first attempt 
to obtain this canted magnetic state from first principles 
in a way that, simultaneously with the atomic potentials 
and effective fields, the directions of the magnetic mo- 
ments are obtained selfconsistently. 

We first performed a calculation for a Pt (111) surface 
in which 8 layers of Pt together with 4 layers of vacuum 
were treated selfconsistently. A seven-atom chain of Co 
together with 10 empty (vacuum) spheres were embedded 
into the topmost Pt layer as schematically indicated in 
Fig. n in order to create a nascent step edge and nested 
Co-chain. Simultaneously, all the nearest neighbors of 
the Co atoms were re-embedded into the respective Pt or 
vacuum layers to allow for relaxation of potentials around 
the Co chain. Thus, an effective embedded cluster of a 
total of 55 atoms was treated selfconsistently. Although 
the number of Co atoms we used in our model chain 
is substantially less than in the experiment, our previ- 
ous experience in calculating magnetic properties of finite 
chains suggests that the local moments and the magnetic 



anisotropy energy contributions of atoms in the interior 
of the chain quickly approach the corresponding values 
of an infinite chain lla . 




FIG. 1: Schematic view of the geometry of a seven-atom 
Co chain along a Pt(lll) step edge. Full circles: Co atoms, 
shaded circles: Pt atoms, open circles: empty spheres. Top: 
side view, bottom: top view of the surface Pt layer with the 
Co chain. The embedded cluster is indicated by solid lines. 
The coordinate systems give reference to the azimuthal and 
polar angles, and 0, that characterize the orientation of the 
magnetization. (Note that in Ref. a different coordinate 
system and the opposite notation for the angles is used.) 

As mentioned before the present implementation of our 
SD scheme serves for searching the magnetic ground state 
of the system. From this point of view it is sufficient to 
consider only the second (damping) term on the right 
hand side of Eq. (QJ. The evolution of the spin orienta- 
tion is then measured on a time scale with a unit (time 
step) of 1/A. A stable ground state is signaled by conver- 
gence of the 9 and (j) angles to a constant and concomitant 
convergence of the constraining fields to zero. In Fig. [21 
the evolution of the 6 and (j) angles is plotted for the first 
100 steps in this artificial time scale for each Co atom in 
the chain. Initially the directions of the atomic magnetic 
moments were set by a random number generator. It can 
be seen that after some oscillations both the 9 and (f> an- 
gles approach a very similar value for all the Co atoms. 
This means that the system rapidly tends to a nearly 
ferromagnetic order due to the strong ferromagnetic ex- 
change coupling between the Co atoms. The initial rapid 
oscillations seen in Fig. [3 are the consequence of the rel- 
atively large constraining fields caused by the large ex- 
change energies that result whenever the individual site 
moments point in very different directions. Once a nearly 
ferromagnetic configuration is achieved, the constraining 
fields drastically decrease, and a rather slow convergence 
of the orientations results. This slow convergence is due 
to the (much smaller) spin-orbit interaction energies and 
it takes about 1000 time steps for the (j) and 9 angles to 
converge. At this point the constraining fields are essen- 
tially zero, as is required for a ground state. The final 
converged state is characterized by a angle of 90°, i.e.. 
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well-known deficiency of the LSDA and often patched 
up using the so-called orbital polarization method or the 
LDA+U method. 
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FIG. 2: Artificial time evolution of the angles 9 (top) and 
the (j) (bottom) angles defining the orientation of the spin 
moments for the seven Co atoms in the finite chain depicted in 
Fig-0 The symbols refer to the following Co atoms numbered 
from the left to the right in the bottom part of Fig. ■ 1, 
O 2, A 3, ♦ 4, A 5, • 6, □ 7. Shown are the resuhs for 
only the first 100 time steps. During the next 900 time steps 
the angles converge very smoothly 



normal to the chain, with a spread of less then 1°, and 6 
angles of 42°. These results are in remarkable agreement 
with experiment |^^, 4> — 90°, 6* — 43°. It is important 
to mention that the obtained ground state is apparently 
not induced by any symmetry of the system. 

Once such a stable state has been established we can 
further analyze the resulting magnetic state in terms of 
the spin and orbital moments. For this purpose their 
size and the corresponding azimuthal angle 9 for each Co 
atom are shown in Table^for each of the Co atoms in the 
chain. The first thing to notice is that the calculated spin 
moments for the inner Co atoms (2.18 fis) are in good 
agreement with the value deduced from experiment (2.12 
P,b) ll^ a nd also with other theoretical studies on infinite 
wires |l7ill9| . Atoms at both ends of the wire have larger 
spin (and orbital) moments than those within the wire 
similarly to our previous findings 0. Although our cal- 
culated orbital moments for the inner atoms (0.19-0.20 
/is) are somewhat larger than the corresponding values 
from other LSDA calculations (0.16 /us jl3 and 0.15 
[l^). they are still much too small when compared to 
the experimental value (0.68 ^b) 0|- This aspect is a 



atom 


Spin moment 


Orbital moment 


moment (/I _g) 


e(deg) 


moment (/is) 0(deg) 


1 


2.23 


41.1 


0.25 39.1 


2 


2.18 


42.5 


0.20 41.5 


3 


2.18 


42.3 


0.19 40.1 


4 


2.18 


42.4 


0.20 41.3 


5 


2.18 


42.3 


0.19 40.2 


6 


2.18 


42.5 


0.20 41.5 


7 


2.23 


41.1 


0.25 39.1 



TABLE I: Calculated magnitudes and orientations of the 
spin and orbital moments in a seven-atom Co chain along 
a Pt(lll) step edge. 

Another interesting feature of the magnetism of finite 
nanostructures is the noncoUinearity of the moments. As 
can be seen from Table U the spin moments of the in- 
ner atoms are quite parallel while those at the end of 
the chain are misaligned by more than 1°. This can be 
traced to larger anisotropy energy contributions at chain 
end sites observed in finite chains earlierfl^. Most in- 
terestingly Table m also reveals differences of as much as 
2° between the orientations of respective spin and or- 
bital moments. This fact underlines the point made by 
Jansen |2Qj, that within DFT the spin and orbital mo- 
ments are aligned only when the ground state refers to a 
high-symmetry direction. 
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FIG. 3: Energy curves calculated using the magnetic force 
theorem for a ferromagnetic seven-atom Co wire at a Pt (111) 
surface step edge as a function of the azimuthal angle, 9. ■ : 
cj) = 90° , A : = 0°. Solid lines serve as a guide for eyes. 

While it is necessary to perform fully selfconsistent cal- 
culations to obtain the details of the non-coUinear ground 
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state, it is interesting to study this result in terms of 
the magnetic force theorem (MFT) [i^j- Assuming fer- 
romagnetic order great simpHfication can be achieved by 
calculating the energy of the system as a function of the 
orientation of a by then uniform magnetization, E(9, 0), 
such that the effective potentials and fields are kept fixed. 
Within this approach only the single particle (band) en- 
ergy has to be taken into account. Fig.Olshows the calcu- 
lated curves of E{0, = 0°) and E(9, (j) = 90°). In these 
calculations we used ground state selfconsistent poten- 
tials and fields obtained from the SD procedure. Clearly, 
the easy axis predicted by the MFT calculations, 9 — 38° 
and (j) = 90° is near the one obtained from the SD calcula- 
tion. Furthermore, the hard axis is obtained a,t 9 = —52° 
and (p = 90°, i.e., by A9 — —90° away of the eas y ax is. 
This is again in good agreement with experiment 

For surfaces and interfaces with uniaxial symmetry the 
easy magnetization axis most often refers to either a nor- 
mal or a parallel direction with respect to the planes. In 
case of a (y, z) mirror plane the energy of a ferromagnetic 
system can be written up to second order as 

E{9,<j)) = En+K2scos29 + K2,2{l-Cos20)cos2(l) 
+ K2,3sm29sm(t) , (4) 

where 7^2, i (« — 1,2,3) are so-called anisotropy con- 
stants. Both calculated trajectories displayed in Fig. O 
coincide almost perfectly with the above function when 
using the parameters K2,i = —0.16 meV, ^^2,2 = —1.06 
meV, and ^^2,3 = —4.81 meV. It should be noted that, 
in terms of these anisotropy parameters, the easy axis 
corresponds to a direction, (j) = '^/'^ and 

'^r^"^"K i^2.r+kJ - 

The anisotropy energy as defined by the energy difference 
between the hard and the easy axes can also be read off 
Fig. 121 the corresponding value of 1.42 meV/Co atom 
compares favorably with the one deduced from the ex- 
periment (2.0 meV) [l^ . It is worthwhile to mention, 
that in terms of similar MFT calculations Shick at al. 
obtained an easy axis along 9 = 18° and an anisotropy 
energy of 4.45 meV 0|. The difference of these val- 
ues with respect to the present results can most possibly 
accounted for in the film geometry and/or the second- 
variational treatment of the spin-orbit coupling used in 
the FLAPW calculations of Ref.:i9]. 

In this study we presented the first application of a rel- 
ativistic ab inito spin dynamics as based on a constrained 
density functional theory to finite magnetic nanostruc- 
tures. In excellent quantitative agreement with the corre- 
sponding experiment, we obtained a canted ground state 
for a Co wire along a Pt(lll) surface step edge. We 
also found that this magnetic state is noncoUinear: a fea- 
ture that is expected to play a key role in nanostructures 



having complex geometry. For the present relatively sim- 
ple geometry the magnetic force theorem proved to be a 
useful tool to interpret the results of spin dynamics cal- 
culations and providing additional information such as 
anisotropy constants and the anisotropy energy. 
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